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55 ■ ABSTRACT 

The determination of the true source polarization given a set of measurements is 
comphcated by the requirement that the polarization always be positive. This positive 
^ ■ bias also hinders construction of upper limits, uncertainties, and confidence regions, 

especially at low signal-to-noise levels. We generate the likelihood function for linear 
polarization measurements and use it to create confidence regions and upper limits. 
, This is accomplished by integrating the likelihood function over the true polarization 

\^ \ (parameter space), rather than the measured polarization (data space). These regions 

' are valid for both low and high signal-to-noise measurements. 



Subject headings: polarization 



Introduction 



^ ' Measurements of linear polarization generally provide values for the degree and angle of polar- 

■ ization. They may also provide upper limits on the degree of polarization even where the nominal 

values P are not much larger than the uncertainties a (i.e. at low signal-to-noise). Polarization 
maps typically plot vectors in regions where measurements are made with high confidence and open 
circles where the polarization is found to be too low to justify a vector, but where useful upper 
limits can still be estimated (e.g. where P + 2a < 1%; Dotson et al. 2000). 

While the distributions of the polarization degree and angle are not normal (gaussian), it is 
still possible to place estimates and uncertainties and confidence levels on their values. The position 
angle distribution, while non-normal, is symmetric about the measured angle so that estimates and 
uncertainties on its value are fairly straightforward (Naghizadeh-Khouei & Clarke 1993). Estimates 
of the polarization degree are complicated by the fact that it must always be positive, introducing 
a bias into any estimate. For any true polarization degree Pq, we expect on average to measure a 
polarization degree P > Pq. This problem has been thoroughly investigated by Simmons & Stewart 
(1985; hereafter SS85). 
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SS85 have also devised a procedure to place uncertainties on the polarization degree. First, the 
expected probability density of polarization (§2) is integrated over a set of possible measurements. 
However, this only places limits on the measured polarization P, the limits varying with the value 
of the true polarization Pq. To correct for this and thus create confidence limits on the parameter 
Pq, SS85 project their results onto the Po-a^is- 

Here we introduce a slightly different procedure for determining confidence regions on the true 
polarization. Rather than integrate the probability density function over the data-space of the 
measured polarization, we integrate over the parameter-space of the true polarization. In this way 
we directly determine the likelihood (§3) that a measured polarization P was drawn from a source 
of true polarization Pq. 



2. Rice Distribution 

Measuring the degree and angle of polarization is equivalent to determining of the amplitude 
and phase of a sinusoidal signal in the presence of noise. The joint and marginal probability 
distribution functions for the amplitude and phase have been discussed in detail by other authors 
(Rice 1945; Vinokur 1965; Simmons Sz Stewart 1985). For completeness we review their results 
here. 

Consider the simple case of rotating a linear polarizer through some angle 9. The resulting 
signal changes as 

s{e) = Pcos2{e + (i)) + SN, (1) 

where the amplitude P is the polarization degree, cf) a relative phase shift (the polarization angle), 
and Sn gaussian random noise. This form can be decomposed into two orthogonal components 

S{9) = {q + qn) cos 2e-{u + UN) sin 29, (2) 

where 

q = P COS 2(f), (3) 

u = Psm2(f), (4) 

gjv is the noise component in phase with cos 26 and uat is the noise component in phase with sin 29. 
In optics q and u are known as Stokes parameters. In electronics they are often referred to as the 
"in-phase" and "quadrature-phase" components, respectively.^ 

Since Sn is gaussian, the measured Stokes parameters q and u are normally distributed about 
the true values qo and uo with equal uncertainties a (i.e. (qj^) = (n^) = a). The values of qo and 



^The optical polarization vector has no direction and is therefore periodic in tt, as is clear from the factor of 2 in 
the arguments to the sine and cosine functions. In electronics the signal is periodic in 2-w by definition, so the factor 
of 2 is replaced by unity in equations (1) - (4). 
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uo are given by equations (3) and (4), replacing P by Pq and (f) by 0o- The uncertainty in the 
polarization P is also given by a. 

The joint probability distribution for q and u is given by the product of their individual normal 
distributions. Transforming to the polar coordinates P and yields 

P2 + p2 _ 2pPq cos 2{(f> - (PoY 



f{P, 



1 

27rcr2 



exp 



2(72 



(5) 



The probability of measuring a polarization in the range [P, P + dP] and [(?!), (j) + d4>] is given by 
f{P, <p) PdPd{^ — <l>o)- Integrating over d{(l) — (f)o) yields the Rice distribution for the polarization 
degree: 

tt/2 



F{P\Po)dP = ^PdP I /(P,</.)d(</>-</>o) 



P 

— exp 



-tt/2 

2(t2 



PPa 



dP_ 



(6) 



where /o() is the zeroth order modified Bessel function. (The pre- factor of | in the integral of eq. 
[6] corrects for the vr-periodicity of (p.) 

The asymmetry of this distribution with respect to P and Pq (evident in Fig. 1 for low values 
of Pq/ct) results in the positive bias of the measured polarization. For high signal-to-noisc the Rice 
distribution is near-normal with mean approaching Pq and a standard deviation approaching a 
(e.g. Rice 1945). 



3. The Likelihood Function 

3.1. Constructing the Polarization Likelihood 

If the true polarization Pq is known, the probability that a measurement will be made in the 
range [P, P + dP], is given by equation (6). However, one is more likely to have knowledge of the 
measurement P and wish to know the probability that the true polarization is within the range [Pq, 
Pq + dPo]. We require the likelihood function L(Po|P) such that the aforementioned probability 
is given by L(Po|P)dPo. For a model characterized by parameters y = {yi,y2, ■ ■ ■ ,ym}, and a 
measured dataset x = {xi,X2, ■ ■ ■ ,Xn}, the likelihood function is constructed from 

n 

L{y\x) = n f{x^\y) (7) 
1=1 

(Eadie et al. 1971). For a single polarization measurement {n = 1) the likelihood function becomes 

L(Po|P) dPo = -^P(P|Po)dPo (8) 
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where the normalization constant is 



= P 

The Hkelihood function is plotted in Figure 2. 



poo 

Nl = / F(P|Po)dPo 
Jo 

fr^'^ (-5) (5) ^ 



3.2. Maximum Likelihood 

A "best fit" estimate of the true polarization can be extracted from the measured polarization 
by utilizing the maximum likelihood. Taking the derivative of L[Pq\P) we have 



dL 



dPo 




^^o(^)=0 (10) 



where /i() is the first order modified Bessel function and P is the maximum likelihood estimate of 
Pq. Complete solutions to this equation can be found numerically (Fig. 3a). At the limits of low 
and high signal-to-noise the solution is 

P = for P/a < V2, and (11) 



P fa y/P^ - (j2 for P/a > 3. (12) 

SS85 have shown that the maximum likelihood estimator does not completely correct for the positive 
bias of polarization measurements. However, of the estimators investigated by SS85, the maximum 
likelihood is superior for low signal-to-noise {Pq/ct ^ 0.7). For higher values, the most probable 
estimator (that which maximizes the F{P\Pq) in Fig. 1) is best (see also Wardle &: Kronberg 1974). 
For Pq/g >4:^ both estimators approach that of equation (12) (SS85). 



3.3. Confidence Regions 

SS85 create confidence regions by integrating the probability density function (eq. [6]) over 
the measured polarization P. To allow for the positive bias in the probability distribution they 
construct confidence regions on Pq by projecting the results onto the Po-axis. However, it is not 
clear that this approach completely accounts for the asymmetry in the Rice distribution. 

To determine the probability that a single measurement P/cr is drawn from a range of true 
polarizations Pq/ct G [Pj/cr, Pu/cr] wc integrate the likelihood function over the parameter space of 
Pq/ct. Following the approach of SS85, limits of measure f3 x 100% for each value of P/ct can be 
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Measured Polarization P/a 



Fig. 1. — Probability density distribution for the measured polarization P/cr, given by equation 
(6). Distributions are shown for several different values of the true polarization Pq/g. 




True Polarization P(,/a 



Fig. 2. — Normalized likelihood function L{Pq\P) given by equation (8). Distributions are shown 
for several different values of the measured polarization P/a. 
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constructed from 



L{P^^\P)A{P^yla) = P + X and 



(13) 








(14) 



Jo 

Choosing A such that (P^ — Pi) is minimized impUes that L{Pi) = L(P„) (for all A > 0) but does 
not generate symmetric confidence regions, especially for low P/a. Prom Figure 3a we see that: 

(1) For measured polarization values P/a < 1.7, 2.4, and 3.0, the lower end of the confidence 

region is bounded by zero for the 68.0%, 95.0%, and 99.0% confidence levels, respectively. It is not 
possible to construct a truly bounded symmetric confidence region for these polarizations; it is only 
possible to construct upper limits (A = 0). 

(2) For polarizations P/a > 4 the confidence regions approach those given by a normal gaussian 
distribution centered on the de-biased value P « ^/W^^ (68% ^ la, 95% ^ 2a, 99.0% -> 2.6a). 

Some authors have placed upper limits of "high confidence" on polarization measurements 
using the criterion that P + 2a < P„ (e.g. Dotson et al. 2000). While this criterion matches no 
single confidence level, it does fall between the 95% and 99% confidence levels even for low P/a 
(Fig. 3b). For determining upper limits at high confidence this criterion seems a reasonable (and 
perhaps more feasible) alternative to solving equations (13) and (14). 

Figure 3b illustrates the use of our results to find confidence regions. To facilitate comparison 
of our results with those of SS85, consider a polarization measurement of P = 2%, a = 0.8%. 
Upper and lower confidence limits for P/a = 2/0.8 = 2.5 are given in Table 1. The differences 
between the two sets of results increase with increasing confidence level. These differences arise 
because the integration over data space and subsequent projection onto parameter space performed 
by SS85 is not equivalent to our integration over parameter space. However, these two methods 
become congruent as the signal-to-noise increases and the Rice distribution approaches a normal 
distribution, [tb] 



The probability distribution of polarization measurements is asymmetric, resulting in a positive 
bias and complicating the generation of uncertainties and confidence levels. We generate confidence 
intervals on the true polarization by integrating the likelihood function over the parameter-space of 
the true polarization, as opposed to the data-space of the measured polarization. These confidence 
levels are valid for both high and low signal-to-noise measurements and approach those given by a 
normal distribution for high signal-to-noise. 



4. Summary 



4 8 4 8 

Meosured Polarization P/a Measured Polarization P/ct 



Fig. 3. — Polarization confidence limits. Grayscale indicates confidence regions of 68.0%, 95.0%, 
and 99.0%. a) Solid curves labeled "upper limits only" and dashed line extrapolations indicate 
regions where the lower polarization limit Pi = 0. At high signal-to-noise the confidence regions are 
bounded by those expected from a normal gaussian distribution (shown as labeled solid lines and 
dashed extrapolations). The dotted white curve indicates the maximum likelihood solution (eq. 
[10]). b) The solid line marks an upper limit which satisfies the condition = P + 2a. Dotted 
lines illustrate the construction of confidence regions for a measured polarization of P/a = 2.5. 



Table 1. Pq/c^ confidence regions for P/a = 2.5 



Confidence Level Simmons fc Stewart 1985 This Work 



(percent) 


Lower 


Upper 


Lower 


Upper 


68.0* 


1.4 


3.3 


1.2 


3.3 


95.0 


0.28 


4.32 


0.02 


3.97 


99.0 





4.96 





4.65 



^Simmons 8z Stewart (1985) use a 67% confidence level rather 
than 68%. 
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